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Abstract 

The compatibility of producing the observed amount of dark matter, as indicated by 
the WMAP data, through the relic abundance of neutralinos with Yukawa unification 
and with the measured rate of b — > 57 is analyzed in mSUGRA and extended SUGRA 
unified models with the inclusion of CP phases. The CP phases affect the analysis in 
several ways, e.g., through the threshold corrections to the b-quark mass, via their effects 
on the neutralino relic density and through the SUSY contribution to the BR(6 — > 57) 
which is sensitive to the CP phases. We present some specific models with large SUSY 
phases, which can accommodate the fermion electric dipole moment constraints and give a 
neutralino relic density in agreement with observations as well as with the b-r unification 
constraint. The possibility of achieving WMAP relic density constraints with full Yukawa 
unification is also explored. 



* Permanent address of T.I. 



1 Introduction 



The Wilkinson Microwave Anisotropy Probe (WMAP) has placed stringent bounds on 
the amount of cold dark matter (CDM) in the universe. The amount of CDM deduced 
from WMAP data is given by|TJ |2] 

tt C DMh 2 = 0.1126™, (1) 

where flcDM = Pcdm/Pc, an d where Pcdm is the matter density of cold dark matter and 
p c is the critical matter density needed to close the universe, and h is the Hubble parameter 
measured in units of lOOkm/s/Mpc. It is reasonable to assume that similar amounts of 
dark matter exist in our Milky Way and in the terrestrial neighborhood, and there are 
many ongoing experiments for the detection of such dark matter in the laboratory. On 
the theoretical side the WMAP data on cold dark matter puts stringent constraints on 
unified models of fundamental interactions since such models are called upon to predict or 
at least accommodate the WMAP data on CDM. As is well known, supergravity unified 
models [3] with R-parity conservation allow for the possibility that the lightest neutralino 
may be the lightest supersymmetric particle (LSP) which could serve as a dark matter 
candidate 2 . A hallmark of many unified models is Yukawa unification. In this paper we 
carry out a detailed investigation of the possibility of accommodating the WMAP cold 
dark matter data in the neutralino LSP scenario but under the constraints of Yukawa 
unification and including the effect of CP phases 3 . Another important constraint is the 
FCNC constraint given by b — > s + 7 which is discussed in some detail in this paper and 
included in the analysis. 

Since the main focus of the analysis is the Yukawa unification constraint on dark matter 
in SUGRA models 4 , we briefly discuss some broad features of this constraint with details 
to follow later. In the supersymmetric framework the unification of the Yukawa couplings 
of the third generation, as predicted in several grand unification models, is rather sensitive 
to the parameters of the SUSY models. Thus, the compatibility of b — r unification at the 

2 There is also revived interest in the possibility that the LSP in SUGRA models could be the gravitino. 
For an update see Ref.[5] 

3 For an analysis of dark matter with CP phases but without inclusion of the Yukawa unification 
constraints see Refs.jBl Q. An analysis of dark matter with quasi Yukawa unification was given in 
Refe-IHIEJ. 

For recent works on dark matter analyses in SUGRA models see Ref.|10| 
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grand unification scale with the observed b and r masses depends sensitively on the sign 
of /i 5 (where \x is the Higgs mixing parameter) as well as on the details of the sparticle 
spectrum [T2*l I13j. Moreover, for most of the available parameter space b-r unification is 
in conflict with other experimental constraints such as the FCNC process b — > 57. The 
more stringent b — t — r unification is predicted in the minimal SO (10) models where the 
quarks and leptons, residing in the 16-plet spinor representation of SO (10), gain masses 
via coupling with a 10-plet tensor representation of S'O(IO) 6 Finally, we 

mention that there can be GUT scale threshold corrections to the Yukawa unification. 
However, typically they are expected to be small [TSj . 

Let us now be more specific and review the situation of Yukawa coupling unification 
in the mSUGRA case with no phases. With universal Yukawa couplings at the grand 
unification scale, the masses of the bottom and the top quark are naturally higher than 
the r lepton mass. This phenomenon arises because of the color interactions which causes 
the Yukawa couplings of the quarks to increase as one goes down to lower energy scales. 
Thus, the running quark masses end up larger than the running charged lepton masses. 
To convert the running mass to the pole mass, one needs to include the supersymmetric 
as well as the standard model (SM) threshold corrections. In particular, it is well known 
that the supersymmetric threshold correction to the bottom quark mass, Am&, can be 
very large. The value of Amj is enhanced for large values of tan f3, where tan (3 is the 
ratio (H u ) / (Hd) and where H u gives mass to the up quark and Hd gives mass to the 
down quark and the lepton. In the mSUGRA case with no CP phases, Amj takes the 
sign of n ^H] (unless the trilinear terms are very large). A negative SUSY threshold 
correction to m& is required in models with b — r unification, in order to obtain a b-quark 
mass in the allowed range. Therefore, b-r unification points toward a negative value 
of the //-parameter. However, a negative fi parameter makes the SUSY contribution to 
BR(6 — > 57) positive and hence it adds to the SM contribution and the charged Higgs 
contribution. As a result, a heavy spectrum is required in order not to exceed the upper 
bound for this branching ratio. 

The SUSY contribution to the muon anomalous magnetic moment also takes the sign 
of fi in mSUGRA |19j . and more generally this contribution depends on CP phases |20j . 
5 We use the sign convention on /1 as in Ref.|llj. 

6 More Higgs multiplets are needed to break the gauge symmetry correctly down to the standard model 
gauge symmetry, but typically these additional Higgs fields do not have couplings to quarks and leptons. 
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However, experimentally the situation is less clear regarding the implications of the — 2 
data. Thus, while the BNL experiment has significantly improved the accuracy of the 

— 2 measurement [21], ambiguities in the hadronic error, which is needed to compute 
the deviation of the observed value from the Standard Model prediction, still persist. Cur- 
rently, the largest source of error in the computation of the Standard Model prediction is 
the 0(a 2 ) hadronic vacuum polarization correction. The most recent evaluation of this 
correction are done by (i) Davier et.al. [22] using the r decay data, and by (ii) Hagiwara 
et.al.j22] using the low energy data from e + e~ — > hadrons. Assuming that the entire 
difference Aa M (where a M is defined so that the effective operator is a fl (e/2m^)fia a /3fiF af3 ) 
between experiment and theory comes from supersymmetry, one finds that the supersym- 
metric contribution for the case of Davier et.al. is (7.6 ± 9.0) x 10~ 10 , while for the case 
of Hagiwara et.al., the value is (23.9 ± 10.0) x 10~ 10 . In this analysis we adopt solution 
(i). In fact, in most of the parameter space we explore the sparticle spectrum is rather 
heavy, and the SUSY contribution Aa^ small, and thus the a M prediction is essentially 
the same as the Standard Model prediction which is consistent with the current data. 
Solution (ii) puts a more stringent constraint. However, the constraint can be softened if 
the universality condition on the soft terms is removed [2U 123 12H| • 

As indicated earlier this paper is devoted mostly to an analysis of the WMAP data 
with the b—r unification constraint. However, we will also briefly discuss b-r-t unification. 
As is well known such a unification requires large tan (3 and for this reason much of the 
parameter space is excluded since it does not correctly break the electro-weak symmetry. 
Several studies has been done, and models such as, e.g., the D-term splitting in SO(10) or 
non- universal Higgs masses can indeed give rise to a viable t-b-r unification |27 | l28 ] l29 ] l30 ] 
|3~T] . However, typically these solutions require a very heavy SUSY spectrum. Thus, the 
predicted dark matter abundance of neutralinos, in models with R-parity conservation, 
will be too high and thus will over-close the universe. As mentioned earlier, models with 
quasi unification have also been investigated jHJ E] • 

In the present work we first analyze the relic density within mSUGRA and show that 
there exist regions of the parameter space where the WMAP relic density constraint, the 
Yukawa unification constraint, and the BR(6 — > 57) constraint can all be simultaneously 
satisfied. We then extend the mSUGRA parameter space retaining universality on the 
magnitude of the soft parameters but allowing non- universality for the phases in some 
sectors. The SUSY contribution to Ami, is phase-dependent [H21 IE] and this allows one to 
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determine the phases in some cases such as to obtain m^Mz) in the experimental range 
and thus achieve b-r unification. Indeed, one finds that with the inclusion of phases b — r 
unification is achievable in a large area of the parameter-space. Further, it is possible 
to find arrangement of phases such that the prediction of the electric dipole moments 
(EDMs) is in agreement with the experimental bounds. The case of full Yukawa unification 
is, however, still (almost) incompatible with the experimental value for BR(6 — > 57) . 
However, it is worth keeping in mind that small flavor mixings in the sfermion mass 
matrices can substantially change the predictions for BR(6 — > 57) while leaving other 
predictions essentially untouched. Thus, in principle, non-zero CP-phases could also allow 
viable b-r-t unification modulo mixings in the squark flavor sector. However, we do not 
pursue this line of investigation in the work here. 

The outline of the rest of the paper is as follows: In Sec. (2) we give a discussion of 
the parameter space of the model and the details of the procedure of the calculations. In 
Sec. (3) we discuss the calculation of the BR(6 — > S7) and resolve some of the ambiguities 
present in the literature in the large tan/3 enhanced contributions, by carrying out an 
independent analysis of the parameters e b (t), e t (s), e&&, which codify these contributions. 
In Sec. (4) we carry out an analysis of the relic density with the b — r unification constraint, 
and with the BR(6 — ► 57) constraint both within mSUGRA and in extended SUGRA 
models with phases. In Sec. (5) we give an analysis of the relic density for the case of the 
full Yukawa unification. In Sec. (6) we discuss the consistency of the analysis with large 
phases with the EDM constraints and give examples of models with large phases consistent 
with the WMAP data, b-r unification and with the EDM constraints. Conclusions are 
given in Sec. (7). 

2 Constraints on SUGRA models with CP-phases 

Within mSUGRA there are only two physical phases, which can be chosen as 8^, the 
phase of the Higgs mixing parameter fi, and «o the phase of the universal trilinear term 
Aq. These phases are severely constrained by the non-observation of the electric dipole 
moments (EDM). The present upper bounds for the EDM of the electron, of the neutron, 
and of the mercury 199 Hg atom are jH3 HIH IHH] 

|4| < 4.23 x KT 27 e cm, \d n \ < 6.5 x KT 26 e cm, C Hg < 3.0 x KT 26 cm , (2) 
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where Cn g is defined as in Ref. |36|. Large phases can be accommodated in several scenar- 
ios such as models with heavy sfermions [37] , models with the cancellation mechanism [SH] , 
models with phases only in the third generation [3H], or models with a non-trivial soft 
flavor structure |3Q]. Here, we use the cancellation mechanism [HE] 7 which becomes possi- 
ble if the SUGRA parameter space is extended to allow for different gaugino phases. The 
model we consider is thus described by the following parameters 



where, & is the phase of the gaugino mass Mj, i = 1,2,3. The value of \n\ is determined 
by imposing electroweak symmetry breaking (EWSB). 

In the analysis we use a top-down approach, and thus impose Yukawa unification at 
the GUT scale, Mgut- For b — r unification we have two independent Yukawa couplings 
at the grand unification scale, i.e., one common /z uni for the b and the r, and one for 
the top-quark. We use these to fit the experimental value of the r and the top masses. 
Unless another value is specified, we fix the top mass at 178 GeV, which is its current 
experimental central value [12] • The value of a s is fixed to be 0.1185. For the r mass 
at the electroweak scale Mz, we use 1.7463 GeV, which takes into account the Standard 
Model radiative correction. Naturally, we also take into account the SUSY correction, 
as derived in [22], when calculating m T . In the case of the full Yukawa unification we 
impose hf, = h t = h T = h uni at M G ut- Therefore, the value of tan (5 is fixed, since the two 
parameters h uni and tan (3 are varied so as to obtain agreement with experimental values 
of m T and m to p- As the b-quark couples to the same Higgs doublet (Hd) as the r lepton, 
its mass is fixed by h un i. Therefore, rrib(M z ) is a prediction of our model and we require 
its value to be within the 2a range, 



as described in [S]. In addition to the above, the other important constraints of the 
analysis are the relic density and the BR(6 — > sj) constraint (see section [3]). 

The procedure for the calculation of the particle and sparticle masses is as follows; 
After choosing a given set of the parameters in Eq.©, we run the renormalization group 
equations (RGEs) down to the SUSY scale, defined as the average of the two stop masses. 

7 For a more complete list of references and for a discussion of the effects of CP phases on low energy 
processes see Ref.|^T]. 



mo, mi/2, tan/5, | A \ : M , a Q , £i, £ 2 , £ 3 , 



(3) 



2.69 GeV < m b (M z ) < 3.10 GeV 



(4) 
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At the SUSY scale the scalar potential is minimized and \fj,\ is calculated along with 
the SUSY threshold corrections to, e.g., the b quark and the r lepton masses and the 
couplings are corrected accordingly. Hereafter, the sparticles are decoupled and the SM 
RGEs are used to run down to Mz- At the electro- weak scale we check if the gauge 
couplings, the Weinberg angle, the top quark and the r lepton masses are in agreement 
with their experimental values. If not, the RGEs are run iteratively until convergence is 
achieved. In the analysis we use the two-loop SUSY renormalization group equations [13] 
except for the trilinear terms, the gaugino and sfermion masses, which are calculated at 
the one-loop level. The SUSY renormalization group equation will also be influenced by 
the CP-phases. However, it is easy to see that neither the phase of the /x-term nor the 
phases of the gaugino masses will run. But, the phases of the trilinear terms run, and in 
general there will be three different phases at the low energy scale namely, a t , and a T . 
at is important as it affects Amj as well as BR(6 — > S7) . However, its value is almost 
fixed by the gluino phase. As shown in Ref.jH], the approximate relation A top oc — M3 
holds at low energy. 

The regions of the mSUGRA parameter space that allows for acceptable relic abun- 
dance can be classified as: (i) the x ~ ? coannihilation region, (ii) the resonance region, 
and (iii) the Hyperbolic Branch/Focus Point (HB/FP) region^!- In a previous work jH], 
we pointed out the strong variation of Am), with CP phases. In that work we focussed on 
the effects induced by the SUSY corrections on the spectrum and their consequences for 
the neutralino relic density. It was shown that the CP-phases have a very large impact on 
the value of the CP-odd Higgs mass Ma, which in turn affects the predicted dark matter 
abundance in the so-called resonance region. The analysis of Ref . jHj used a bottom-up ap- 
proach by fixing the value of m^iMz) to its central value. In this work we use a top-down 
approach and large effects of the CP phases are not seen. In fact, the predicted neutralino 
relic abundance, turns out almost independent of the phases in the resonance region. In 
the stau coannihilation region there is also very little dependence on the CP-phases, ex- 
cept for the trilinear phase. As we show below, the HB/FP region cannot give rise to 
Yukawa unification within our model. In the calculation of the relic density we take into 
account the CP even-CP odd Higgs mixing. In the MSSM, after spontaneous breaking of 
the electro-weak symmetry one has at the tree-level two CP even Higgs (h°, H°) and one 
CP odd Higgs (A). In the presence of CP violating phases these mix, producing mass 
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eigenstates (H®, H®, H®), which are no longer eigen-functions of CP [dT] 8 

The most important supersymmetric threshold correction is the one to the bottom 
mass. At the loop level the effective b quark coupling with the Higgs is given by 



- C bbH o = (h b + 5h b )b R b L H? + Ah b b R b L H? + H.c. (5) 
The correction to the b quark mass is then given directly in terms of Ah b and 5h b by 

Am b = [ite(^) tan/3 + (6) 
h b h b 



We use the full analysis of Am b derived in [32] • The largest contributions to Am b are 
the gluino and the chargino exchange contributions. The gluino exchange contribution is 
proportional to M^fi, and will therefore depend on the phase combination 9^ + £3. The 
chargino exchange contribution is usually smaller, except for very large values of \A t \, 
since it is proportional to A t \x. Its dominant phase dependence is given by 6^ + a t , and 
it has the opposite sign of the gluino contribution in a large region of the parameter 
space. When evaluating h b at M SU sy, we take into account threshold corrections using 
the relation 9 

The SM Yukawa coupling is evolved down to the electroweak scale, and the bottom quark 
mass 

m b (M z ) = h S b M ^ cos (3 , (8) 

is calculated and compared with experiment. Similar expressions hold for the r lepton 
with b replaced by r. For the top quark at the Z scale one has 

m t (M z ) = -^sin/3/i t SUSY (l + Am t ) (9) 

where 

Am t = [Re(^) cot/3 + Re A]. (10) 
fit <H 

A full analysis of Am t is given in Ref. jH2j- However, in the region of interest which 
corresponds to large tan /3 the correction to the top quark Yukawa is essentially negligible. 

8 For further details regarding the implications of these CP even-CP odd Higgs mixings on neutralino 
dark matter analysis see Ref.|Hj. 

9 This relation resums the SUSY self-energy leading order logarithmic corrections |47j . 
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3 BR(6 -> S7) with CP-phases 

The present average for the BR(6 — > 57) derived from the available experimental data jlHj 
is found to be, 

BR(b - 87) = (3.54j£5j) x 10- 4 , (11) 

by the Heavy Flavor Averaging Group [IB]- The error includes an uncertainty due to the 
decay spectrum as well as the statistical error. The theoretical SM prediction is |5T)| I5T] , 

BR(b -»• S7) = (3.70 ±0.30) x 10" 4 . (12) 



The above result uses the MS running charm mass instead of the pole mass. It was 
claimed in Ref. [SU] that this consideration reduces the NNLO uncertainty in the SM 
calculation. However, other analyses [221 EH] question the theoretical precision of Eq.(|12|) 
predicting a lower central value for the SM result. In any case, the result of Eq. (jl2|) 
appears to be a good benchmark value for the SM prediction to work with. 

The dominant SUSY contributions from the charged Higgs exchange include the tan f3 
enhanced NLO corrections, which contribute to the Wilson coefficients CV and C% (these 
are coefficients of the operators 7 = j^m b (s L o-^ u b R )F^ v and 8 = -^m b (s L o- t , u T a b R )G c [ l 
These contributions can be codified in e' b (t), e t (b) and e bb which enter in the Lagrangian for 
effective interaction involving the charged Goldstone boson and the charged Higgs boson 
as follows 

c = g^^-S^ tt^Tp " A} (13) 

+ -7^H + {T + £ m.V.At. . . } + H.c. 

V2M W ^ tan/3 V l + e£ 6 tan/3 

where Vi 3 - is the CKM mixing matrix. Evaluation of e' b (t), e t (b), and e bb exist in the litera- 
ture j33JESl; but there is some ambiguity concerning the signs of some of the terms among 
the above groups. To resolve this we carry out an independent analysis of these quantities 
for the same loop diagrams as in the previous works, including also their dependence on 
CP phases, which was taken into account only in one analysis previously. Our analysis is 
derived using the work of Ref. |56j. We find 



2 2 n 

<£(*) = "EE ^^Dl 2j D m [— cot j3A t D bXj D* t2i + ^D b2j D* tll + m t cot /3D b2j D* m 

i=l j=l 671 m b 



8 



m 2 t m 2 w . „ „, m ? m V 



cot (3D blj D* tli ^ sin /? cos /3D blj D* tli ] ,#( 

7716 7716 |77lg| l m 9l 

4 2 2 

+2EEE1- c °t PAD blj D; 2i + /xDwyD^ + m t cot /3A,2; A* 2 * 

fc=l i=l j=l m?l 
2 2 

Tn Tn 

cot (3D hlj D* tli ^ sin cos pD blj D* tll ] (14) 

777-6 ?77f> 

,(al k D h - 7iO; 2i )(&A 1( + a^J^-i-ir^, ^) 

In the above D q is the matrix that diagonalizes the squark mass 2 matrix M|, i.e., 

Dt M |D g = diag(Mj,M?) (15) 

and H(a, b) is defined by 

H(a,b) = - £ 7rlna + - — -ln& (16) 

V ; (l-a)(a-6) (1 -&)(&- a) V ; 

where a b k, Pbk,1bk for the b quark and the corresponding coefficients for the t quark are 
as defined in Ref-jHIj- Similarly for e t (s) we find 

e t (s) = EE ;i! /;; ; ,/V— tan/M^L^. + v*D sXi D* t2j + m s tan /3L» s2iJ D* 2 . 
i=i j=i m * 

777^ 777^ 1 777^ 777 ~ 

+-*taa/3£> 8li £>; y - sin /3 cos f3D sli D* ] ] H{-^-,-\) 

m t J m t J \m§\ \fng\ \ m g\ 

4 2 2 

771 

"2 E E Et— tsmPA*D s2i D* tlj + ^D sU D* t2j + m s tan /3D a2i D^ 

k=X i=l j=l m * 

2 2 

+^ta,n(3D sU D* - ^ sin (3 cos (3D sli D*} (17) 
m 4 J m t J 

1 1 „,< 



x(P*sk D *sii + a«fc^* 2i )(a tfc D t y - 7t*Asy) jg^ZTT #( 



2 ' ™2 



777 ^0 777 o ?T7 



Finally, our analysis of e b b gives 

A A 2a s _i 6n * rM 2 m w cos/3 1 1 . 2 . 

e w> = -LL^ e ^h- ^ T~ i -o + o sm Ow)D b uD blj 

37T m fe cosfcV 2 3 J 



In 1 ™l 

--sin 2 0^D 62i i^} sin/3 + //Ai^l^^H^ 
3 |m§| |^g| l m sl 
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V^V^V" 2t M zmw cos 1 2. 2 2 . 2 . 

2^ 2^ 2^ [ {(q ~ q sm e w)D tli D tlj + - sin e w D t2 iD t2j } sm(3 

7T}/fo COS t/p^/ 



=i j=i fe=i 



m 6 m fe 



^ cot (3{D m D* tlj + D t2i D* 2j } - — cot (3A* t D t2 iD* tlj \ 



,2 /m? 



The form factor i?(a, 6) in the above equation can have a = b and in this case it reads 

H(a,a)= . V.J l-a + lna] (20) 
(a - 1)2 

Before proceeding further we give a brief comparison of these results with the results of 
the previous works. The analysis of e 6 (t) may be compared to e#, of Ref.jHZj in the limit 
of large tan /3 and small squark mixings. In this case the limit of the first two lines in 
Eq.(|15|) agrees with the result of Ref. [27]. However, the limit of the last three lines of 
Eq. (fTK|) have an opposite sign to that of Ref . jHZj ■ Here our analysis is in agreement with 
the result of Ref. [SB]- 

Next we give a computation of e t (s). Approximating Eq. ffTHj) we find 



i= i i= i 3tt |m g | \m g \ 2 \m g \ 2 

h 2 s A* 22 m l 
~ : „ s 9 — —X^kXik\D S 2i\ \D t ij\ H(- % — -,- i--) (21) 

The analysis of Ref. [HZ] computed only the first line of Eq. (|21|) and for this part we agree 
with their work when we take the large tan f3 limit and the limit of small mixing angles 
of our result. The work Ref. [HI] gives results corresponding to Eq. (j21|) . However, here we 
find that we have a disagreement with the sign of the second part of their Eq.(16). 

Our analysis of given by our Eq.(fT§]) agrees with the analysis of Ref.jHZI in the 
limit of large tan (3 and in the limit of small squark mixings and here there is a general 
agreement (taking account of typo corrections) among various groups in the limit of no 
CP phases. Our analysis like that of Ref. jHZI takes into account the full dependence on CP 
phases. In the numerical analysis to be presented below, we have used the code provided 
by micrOMEGAs jHH] in the CP conserving case. This code agrees with the codes used 
by other groups jSH] • In the CP violating case we have combined the codes of Refs. [58] 160] 
with our own codes of the SUSY contributions with CP phases. 
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For the uncertainty in BR(6 — > sj) we use a linear combination of the errors on 
Eqs. (HH) and (fT2Jl. At the 2-a level one has 

2.3 x 10" 4 < BR(b -> S7) < 4.7 x 10" 4 . (22) 

The numerical analysis given below is controlled essentially by the upper bound in Eq.([22|). 
In order to obtain the correct value of rrib(mz) one needs the phase combination 9 ^ + £3 
to be close to 7r 10 . In this case the chargino contribution to the BR(6 — > 37) is positive 
and therefore the lower bound is not reached for the values of the SUSY parameters in 
our study. We turn now to the details of the numerical analysis. 

4 WMAP Dark Matter Constraints and b — r Yukawa 
unification 

It is useful to first summarize our results in the mSUGRA case, where all CP phases are 
either zero or ir. The relation = h T can be satisfied for a wide range of soft masses 
in the MSSM with real universal soft terms. To discuss the dependence of m^iMz) on 
tan/3 we consider two representative set of soft parameters: (i) mi = 800 GeV, A = 0, 
m = 300 GeV, and (ii) mi = 800 GeV, A = 0, m = 600 GeV. In Fig. (JTJ) we study the 
H > and \x < cases for each set. The lines corresponding to case (i) are interrupted 
when nif < m x , while for case (ii) large values of tan/5 are incompatible with EWSB. 
We include a reference line ignoring the SUSY threshold corrections (i.e., Awij = 0). 
Fig. (Q) exhibits the well known phenomenon, that Awij is positive for fi positive and 
therefore the theoretical prediction for the b quark pole mass is too high, lying outside 
the experimental range. Thus b-r unification does not occur in this case. When fi < 0, 
on the other hand, Ami, is negative and the theoretical prediction for the b quark mass 
can lie within the experimental range for values of tan/3 between roughly 25 and 45. A 
similar analysis of rrib(Mz) but as a function of m is given in Fig.Q. Here we consider 
only the fi < case and find that the theoretical prediction of nibiMz) can lie within 
the corridor allowed by experiment for a range of m values. However, one finds that for 
very high values of mo, i.e., for values above 5 TeV and beyond, a region which includes 

10 The phase combination is drawn to smaller values for large tan /3 and for the full Yukawa unification 
it ends up close to 7r/2. 
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the Hyperbolic Branch (HB) /Focus Point (FP) region, Yukawa unification is not achieved 
with universal soft parameters. 

We extend the analysis now to include the relic density constraints. Fig.© shows 
an area plot in the m — mi plane and all the three interesting dark matter regions in 
mSUGRA can be seen. In Fig.(jnj) the coannihilation area and the resonance area overlap. 
The HB/FP area, which is the region adjacent to the area with no EWSB (/i 2 < 0), is 
incompatible with any kind of Yukawa unification within the framework of universality of 
soft parameters at the GUT scale as already seen in the analysis of Fig.(J2J). 

The HB/FP region moves to lower values of m as tan(3 and m t decrease and this 
variation, especially with m t , can be very large (we present our analysis in Figs. (j2l3|) 
using m t = 176 GeV in order that the HB/FP area appear below 20 TeV). With m t at its 
lower bound the HB/FP region appears at values of mo of 5-6 TeV, but here m;, is already 
too high. With large values of tan/3, m\ becomes negative before \x becomes small and 
therefore there is no inversion of the gaugino/Higgsino components in the composition 
of x°- Correspondingly, the HB/FB is not reached. This is the case for the line of 
tan/3 = 48 in Fig|2j Therefore, overlapping of the HB/FP region and the allowed m& area 
is not possible. Moreover, phases cannot improve the situation as Awij is very small in 
the HB/FP (below 5%), and thus cannot lower the value of m& sufficiently. 

In Fig. (j3J) we further analyze the mSUGRA case with area plots in the mo — mi 
plane for four values of tan/3: 30, 35, 40, 45. Fig. (JJJ) shows that the relic abundance, 
b — r unification, and b —>■ 57 constraints can be simultaneously satisfied for a narrow 
range of parameters for values of tan/3 in the range 30-45. The BR(6 — > sj) constraint 
is a major restriction, since both SUSY and Higgs contributions to the branching ratio 
add to the one from the Standard Model, and thus one needs a relatively heavy spectrum 
such that the BR(6 — > sj) prediction remains below the experimental upper bound. The 
b-r unification constraint and the WMAP constraint further reduce the parameter space. 
Even so, one finds that there exist regions of the parameter space for all the four cases 
in Fig. (j3J) consistent with the WMAP data under the b — r unification and b — > 57 
constraints. 
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4.1 Effects of CP phases 

To determine the impact of phases on the above picture we choose two representative 
points from Fig. HJ 

a) tanp = 30, m = 290 GeV, M 1/2 = 800 GeV, A = GeV (23) 

b) tanp = 40, m = 710 GeV, M 1/2 = 800 GeV, A = GeV (24) 

These points are chosen because in the absence of phases the WMAP relic density con- 
straints are satisfied by different mechanisms for these two cases. Thus, for the point in 
Eq.()23J) the WMAP constraint is satisfied due to x ~ ? coannihilations. In contrast, for 
the point in Eq. (j24J) the WMAP constraint is satisfied due to a resonance in the Higgs 
mediated annihilation of \ — x- To determine the effect of phases we study the most rel- 
evant phases for the processes that we consider. The phases £i and £2 have little impact 
on b — > S7 and Am&. For simplicity we set them to zero in this section. We have already 
discussed the phase combinations that play an important role in the analysis of Amj. 
For the analysis of b — ► 57 we find that the same phase combinations, i.e, Arg(/j,A t ) and 
Arg(fiMs) are the important ones. 

We now discuss the specifics of the point in Eq. (J23|) . which as already stated is in 
the x ~ t coannihilation region. In Fig.fjSJ) we analyze the BR(6 — > 57) and the b — r 
unification constraints in the 9^ — £3 plane, and as is seen the point satisfies the b — > 57 as 
well as mj (Mz) constraint in the mSUGRA case with a negative /1. The figure illustrates 
that the inclusion of phases changes the value of m^iMz) and BR(6 — > 57) drastically. 
However, the two above mentioned constraints have a tendency to conflict with each 
other, even with the inclusion of CP phases. Nevertheless, we find that there exists a 
substantial overlap of the areas allowed by the bounds on m^iMz) and BR(& — > 57) . 
At the same time the predicted Qc'DAih 2 remains inside the WMAP bounds because the 
phases do not affect significantly the ratio m x /m,f and hence the relic density prediction 
remains dominated by coannihilations. It is also instructive to study the effects of a®. 
In Fig. © we analyze the dependence on \Aq\ and «o for the point £3 = 0.3 rad and 
6* M = 2.4 rad of Fig.©. In the analysis of Fig. (jnj the ratio mf/m x does not exceed 1.08, 
and thus we remain in the coannihilation region allowing for the satisfaction of the relic 
density constraints. Furthermore, it is also possible to satisfy BR(6 — > 37) and mj, bounds 
simultaneously for a wide range of \A \ and a . 
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Next we analyze the implication of phases for the point in Eq. (|24[). As already 
indicated this point is within the resonance region in the mSUGRA case. However, the 
point produces a value for the BR(5 — > S7) outside the experimental bounds as may be 
seen from Fig. (0J). The effect of varying £ 3 and 9^ is analyzed in Fig. (J7|). Here one finds 
a substantial overlap of the areas allowed by the bounds on m b (M z ) and BR(6 — > 57) 
while the relic density prediction remains within the WMAP bounds. As already stated, 
this analysis is substantially different from the one given in Ref. |x at tan/5 = 40. There 
m^iMz) was fixed and the dependence of Amj on the phases has a big effect on the 
resonant channels. For the present case, the bottom Yukawa has only a small fluctuation 
due to the unification condition at the GUT scale. Thus its effect on the Higgs mass 
parameters through the RGE's is not as large as the one found in the analysis of Ref. jH] . 
Thus in the analysis of Ref. jH] no unification condition was assumed for the Yukawa 
couplings, and the only requirement on them was to predict fixed values for the fermion 
masses. In the case of m^, the effects induced by the phases via Am), were compensated 
by variations on hb so as to obtain a fixed rrib(Mz)- Since such adjustments of hb induced 
large changes on the Higgs parameters, the relic density was very sensitive to the phases. 
In the present case, hb is approximately fixed by the condition hb = h T at the GUT scale 
(where h T is determined by tan/3 and m T ). 

The fluctuation of Amb with the phases enters in the prediction of mb{Mz) which is 
allowed to vary in its experimental range. Thus in contrast to the analysis of Ref. [S] hb is 
not adjusted as the phases vary in the present analysis. The value of hb is approximately 
fixed by the condition hb = h T at the GUT scale. Consequently, the phases do not have 
a big effect on the Qh 2 prediction in the present scenario. For example, in Fig. ((ZJ), flh 2 
varies only in the approximate range 0.10 - 0.13. The effects of variations with for 
the point Eq. (J21J) are analyzed in Fig. (JHJ). Specifically, Fig. (jHJ) gives an analysis of 
the neutralino relic density in the \A\q — «o plane for the input of Eq. (ffilj) along with 
9fj, = .5, £3 = 1.7. One finds a considerable structure here exhibiting the important 
effects of «o in this case. The relic density remains within the WMAP bounds, in the 
dark hatched area, while the area above the dashed line has a relic density below the 
lower bound. 
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5 WMAP dark matter constraint and full Yukawa 



unification 

In the above we discussed the satisfaction of the WMAP relic density constraints consis- 
tent with the BR(6 — » s + 7) and b — r unification constraints within mSUGRA and its 
extensions including phases. It was seen that the loop corrections to the b quark mass 
(and to the r lepton mass) play an important role in accomplishing b — r Yukawa unifi- 
cation at the GUT scale consistent with the experimental values for the r lepton and the 
b quark masses. The values of tan j3 used in the analysis above were fairly large, lying in 
the range up to 40 — 45. When tan f3 exceeds these values the possibility that full Yukawa 
unification for the third generation holds becomes feasible. Here we investigate this pos- 
sibility in further detail to determine if WMAP relic density and the BR(6 — > s + 7) 
constraints can also be simultaneously satisfied. In the analysis we will allow for the 
dependence on CP phases. 

In Fig. Q we present an analysis of full Yukawa unification and we also display the 
constraints of relic density and of BR(6 — > s + 7). We impose full Yukawa unification at 
the GUT scale, the value of tanf3 is therefore fixed by the experimental r and top masses. 
As before, m^iMz) is a prediction. Typically, there are two main constraints on m and 
mi/2 for a given A Q . These are the condition of radiative EWSB (or almost equivalently 
rrtA < 120 GeV) and the condition that the LSP be neutral. The constraints on m and 
mi/2 such that both conditions are met were described in an early paper ^3] and also 
emphasized in Ref . |62j which gave the relation 

m\ = otm\i 2 — — constant (25) 

where the coefficients a and (3 are positive and ~ 0.1, and the constant is ~ Mz- Thus for 
fixed rriA one has a hyperbolic branch. Furthermore, the requirement that LSP be neutral, 
i.e. m x < m^, makes another cut in the allowed area. While Fig. Q exhibits a narrow 
area where the WMAP relic density constraint is satisfied, one finds that m^iMz) is 
outside the experimental bounds (the line rrib(Mz) = 2.50 GeV is presented as a reference 
line in Fig. ©). 

In the whole figure the value of rrib{Mz) lies below the lower experimental bound. 
The region satisfying the BR(6 — » 57) bounds is also exhibited in Fig. (JHJ). The charged 
Higgs contribution is enhanced in this case due to the low values of its mass, (lines 
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corresponding to the values = 300 and 500 GeV are given as reference). Since 

the SUSY contribution is also positive the value of BR(6 — > s + 7) lies below its upper 
experimental limit only for the small region found at mi/2 ~ 2900 GeV. For the SM 
contribution we followed the considerations of JSOj by using the MS running charm mass, 
so that ^ = 0.29. In this case the central value of Eq.^Jfor the SM prediction is obtained. 
However, as argued in Ref. |52[ IB^] the theoretical SM prediction is possibly lower. Thus 
as an illustration we also give an analysis using the pole mass ratio ^ = 0.29 which leads 
to a SM prediction of 3.33 x 10" 4 . 

We investigate now the implications of extending the parameter space by CP phases 
for a selected point in the coannihilation region. Fig. (fTUj) shows the 9 ^ — £3 plane for 
m = 880 GeV, M 1 / 2 =1500 GeV and A = 0. The value of tan/3 lies in the range 51 - 54.5. 
The prediction for the relic density remains in the WMAP range, since the neutralino 
remains in the coannihilation area. The regions where rrib(Mz) and BR(6 — » 37) lie 
inside the experimental bounds are shown. There is only a rather tiny region, roughly at 
6* M = 7r/2 and £3 = 0, where the Yukawa unification constraints and the BR(6 — > 57) are 
simultaneously satisfied. This area is affected by the uncertainty in the determination of 
the SM value for BR(6 — > 57). The area is significantly enlarged when the ratio ^ = 0.29 
is used in the BR(fe — > S7) computation. 

6 Consistency with the EDM constraints. 

With inclusion of phases, one has to account for the satisfaction of the EDM constraints. 
In the following we demonstrate that, there exist regions in the parameter space, where 
the WMAP, the b — r unification, the BR(6 — > 57) as well as the EDM constraints are all 
satisfied when the phases are large. In Table 1 we define two points, one for tan/3 = 40 
and another for tan/3 = 45 where all constraints are satisfied as shown in Table 2. 

A more detailed exhibition of the value of m^iMz) and BR(6 — > s + 7) as 9^ and £3 
varies is given in Figs, (jllj) and (|T2"|) while £1, £2 are set at the values given in Table 
1. It was shown in Ref. |63j . that if the EDM constraints are satisfied for a given point, 
there exists a scaling region, where mo, mi/2, Aq scale by a common factor A, in which the 
EDM constraints also are satisfied, for a reasonable range of A around one. The allowed 
range of A depends on other dynamical parameters. For Point (i) m and m 1 / 2 are related 
by m = 0.832 • while for Point (ii) this relation becomes m = 1.80 • m^. For 
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point (ii) the EDM constrains are satisfied down to m\/2 ~ 750 GeV. In Fig. (|13|) we 
display BR(6 — > 57) and the neutralino relic density for the case of two points in Table 1. 
The analysis shows that m^iMz) remains inside its experimental range for the range of 
parameters shown in the figure. The qualitative behavior of the relic density in both cases 
can be understood by comparison with the corresponding cases in Fig. For tan (3 = 40 
the line mo = 0.832 • mx/2 has a sizeable overlap with the WMAP area, whereas for 
tan/? = 45 the line mo = 1.832 ■ mi/2 intersects the WMAP area. The values of m^iMz) 
ranges from 2.80 to 2.86 GeV for the case (i) and from 2.84 to 2.96 GeV for case (ii). 



Point 


m 


mi 


A 


tan/3 




a A 


£1 


6 


£3 


(i) 


1040 


1250 





40 


2.9 





1.0 


0.15 


0.5 


(ii) 


1980 


1100 





45 


0.6 





0.5 


-0.6 


1.6 



Table 1: Values of the parameters for point (i) and point (ii). 



Point 


\d e \ e.cm 


\d n \ e.cm 


C Hg cm 


Qh 2 


BR{b -> 57) 


m b {M z ) 


(i) 


1.33 x 10~ 27 


8.87 x 10~ 27 


1.72 x 10~ 26 


0.099 


4.44 xlO" 4 


2.85 


(ii) 


1.87 x 10~ 28 


2.71 x 10" 26 


1.13 x 10- 26 


0.112 


4.37xl0" 4 


2.92 



Table 2: The EDMs, the relic abundance Qh 2 , the branching ratio BR(6 — > 57) , and the 
mb(Mz) prediction for point (i) and point (ii) as defined in Tabled 

7 Conclusions 

The main focus of this work is an analysis of the neutralino relic density consistent with 
the WMAP data under the constraint of b — r Yukawa unification, and the constraint 
of b — * s + 7 branching ratio. In the analysis of the 6 — * s + 7 branching ratio, we 
have included the tan f3 enhanced NLO corrections which contribute to the Wilson co- 
efficients. These enhancements are codified via the epsilon terms defined in Eq.(13). 
There is ambiguity in the sign of some of the terms among the various groups. To 
resolve this we carried out an independent calculation of these quantities as discussed in 
Sec. 3. The analysis is carried out within SUGRA unified models where universality on 
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the magnitudes of soft parameters at the GUT scale is assumed, but we allow for CP 
violating phases and specifically allow non-universality of the phases in the gaugino mass 
sector. First we give an analysis for the case when all the soft parameters are real. This 
is the mSUGRA case, and here we find that for values of tan f3 in the range 27-48, one 
obtains an amount of dark matter consistent with WMAP as well as consistency with b — t 
unification and with the b — > s + 7 constraint. An interesting phenomenon that arises is 
the following: There are three regions in the mo — mi parameter space where relic density 
and the b — > s+7 constraint can be satisfied in general. These consist of the coannihilation 
region, the resonance region, and the HB/FP region. Of these only the first two can satisfy 
the Yukawa unification constraint. Thus the constraint of Yukawa unification narrows the 
available parameter space by eliminating the HB/FP region. We then extend this analysis 
to include phases and show that new regions of the parameter space allow for consistency 
with the WMAP data and other constraints extending the allowed region of the parameter 
space. In b-r unification case we find explicit phase arrangements such that the EDM 
bounds are satisfied, m^iMz) and the rate for BR(6 — > 57) lie within their experimental 
ranges, and the prediction of the neutralino relic density lies within the WMAP bounds. 
We have also given an analysis of the full b — r — t Yukawa unification constraint with 
inclusion of CP phases. We find a small area where rrib is predicted inside the experimental 
range and the BR(6 — > 57) bound is satisfied. Furthermore, the relic density of neutralinos 
lies within the WMAP bounds due to % — f coannihilations. However, this area is rather 
small and moreover we could not find phase arrangements satisfying the EDM constraints. 
It is conjectured that inclusion of additional non-universalities could rectify the situation. 
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Figure 1: The value of m^iMz) versus tan/3 assuming = h T at the GUT scale with 
m 1/2 = 800 GeV, A = GeV, m = 300 GeV (solid lines), m = 600 GeV (dashed lines). 
The thick (thin) lines have /i < (// > 0). The dot-dashed line is plotted with Am b = 0. 
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Figure 2: The value of m^iMz) using the constraint of b — r Yukawa unification in 
mSUGRA. Furthermore, m t = 176 GeV, /x < {B^ = tt), m 1/2 = 800 GeV, A = and 
the value of tan f3 has been varied as indicated on the curves. Lines ends at values of mo 
where EWSB is no longer satisfied, except for tan (3 = 48, where the CP odd Higgs mass 
is too small. 
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Figure 3: Analysis of neutralino relic density, BR(6 — > 57) and of m^iMz) including the 
HB/FP region for tan/3 = 40, A = and m t = 176 GeV. Areas contoured by the dashed 
line has a neutralino relic density is inside the WMAP bounds. The area above the solid 
line predicts m^iMz) > 3.10 GeV while the area inside the dashed (dot-dashed) line is 
excluded by the lower bound on rrih (the upper bound on BR(b — > 7)). On the lower dark 
shaded area m x > nif while on the upper EWSB is not achieved. The thiner dashed line 
indicates m Y + = 103 GeV. 
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Figure 4: Analysis of the neutralino relic density with the b — r Yukawa unification 
constraints in the mo — M\i<i plane when the soft terms are universal and real with \x < 
(Ofj, = 7r), Ao = 0, m t = 178 GeV,and tanfi = 30,35,40,45. The lines and shaded areas 
are as described in FigJHl The area inside the solid line for the case tan(3 = 45 predicts 
m b (M z ) < 2.69 GeV. 
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Figure 5: A plot in the £3 — M plane for tan/? = 30, m = 290 GeV, mi/2 = 800 GeV, 
Ao = and £1 = £2 = 0. The area inside the solid lines predicts rrib(Mz) within the 2-a 
experimental range. The area inside the dot-dashed line is excluded by the BR(6 — > 57) 
constraint. 
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Figure 6: Analysis of b — r unification in the {Ao/min] — plane with phases corre- 
sponding to the point in Eq.(j23J) with £ 3 = 0.3 rad, 9^ = 2.4 rad and £i = £2 = 0. In 
the dark shaded area, m x > m f while within the ruled area bounded by the dot-dashed 
line BR(6 — > 57) exceeds its upper limit. The grey area predicts m^iMz) inside the 
experimental bounds and fills the whole plane. 
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Figure 7: Same as FigJHl except that tan(3 = 40, m = 710 GeV, mi/ 2 = 800 GeV, 
A = £1 = 6 = 
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Figure 8: Analysis of the neutralino relic density in the [Ao/mi \ —cto plane for £3 = 1.7 rad, 
6* M = 0.5 rad while the other parameters are the same as in Fig[7| The area marking are 
the same as in Fig|BJ The dark hatched area, contoured by the dashed, line predicts Qh 2 
inside the WMAP bounds. 
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Figure 9: Analysis of the neutralino relic density, of BR(b — > s + 7) and of b — t — r 
unification in mSUGRA when \x < (9^ = 7r)and A = 0. On the upper ruled area 
rriA < 120 GeV, while in the lower dark shaded area the lightest neutralino is not the LSP. 
The narrow area bounded by dash lines corresponds to the WMAP favored relic abundance 
prediction, while in the area, bounded by the dot-dashed thick line, the prediction for 
BR(b — > S7) is inside the experimental bounds. The thin dot-dashed line indicates the 
expansion of BR{b — > 37) allowed area when the ratio ^ = 0.29 is used in computation 
of the SM contribution. For the range of parameters exhibited, the prediction of rrib(Mz) 
is below the experimental bound, the solid line corresponds to a prediction of m^Mz) = 
2.50 GeV. The double-dotted-dashed line corresponds to the indicated values of mjj+. 
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Figure 10: The full b — t — r unification allowed areas (bounded by solid lines) and 
BR(6 — > s'-f) excluded areas (bounded by dotted-dashed lines) in the 9^ — £3 plane, for 
m = 880 GeV, Mi/ 2 = 1500 GeV, A = GeV and all the remaining phases are set to 
zero. The area within the thin dotted-dashed line shows the expansion of the b — > S7 
allowed area when — = 0.29. In the dark shaded area the lightest neutralino is not the 
LSP, while in the dark hatched area in the right corner < 120 GeV. 
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Figure 11: Analysis of the b — r unification and the BR(6 — > 57) constraints in the 
6* M — £ 3 plane for tan/3 = 40, £1 = 1.0, £2 = 0.15 and A = 0. m and mi/2 satisfies 
mo = 0.832 • mi/2- In the area contoured by solid lines m& is inside the experimental 
range, while the area inside the dotted-dashed line is excluded by the BR(6 — > 57) bound. 
The thicker (thiner) lines correspond to mi/2 = 1250 GeV (1050 GeV). 
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Figure 12: Analysis of b — r unification and of BR(6 — > 57) in the Q yi — £3 plane for 
tan/3 = 45, £1 = 0.5, £2 = —0.6 and Aq = 0. mo and mi/2 satisfies the equation 
mo = 1.8 ■ mi/2- In the area contoured by solid lines m& is inside the experimental 
bounds, while the area inside the dot-dash line is excluded by the higher(BR(6 — > 57) 
bound. The thicker (thiner) lines correspond to m^j = 1000 GeV (800 GeV). 
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Figure 13: Analysis of the relic density and of BR(6 — > S7) for the points in Table 1. 
For the case (i), m is constrained to satisfy m = 0.832 ■ m 1( / 2 and in the case (ii), 
m = 1.8 • m 1/2 . 
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